function y=beta1(x)
y=x.^3-x.^2;
end